EXISTENCE OF A UNIQUE STRONG SOLUTION TO THE DMPK EQUATION 
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ABSTRACT. For the transmission of electrons in a weakly disordered strip of material 
Dorokhov, Mello, Pereyra and Kumar (DMPK) proposed a diffusion process for the 
transfer matrices. The correspoding transmission eigenvalues satisfy the DMPK sto- 
chastic differential equations, like Dyson Brownian motion in the context of GOE/ GUE 
random matrices. We control the singular repulsion terms of this SDE with a stopping- 
time argument, and its degenerate initial condition by an approximation procedure, and 
thereby establish the DMPK equation to be well posed. 



1. INTRODUCTION 



o 

£ 

h— 5 F. J. Dyson [1 1 J introduced an N x N matrix-valued Brownian motion, to be precise 

00 a stochastic process (A"2(t)) t>0/ 

(1) X 2 (t) = Y 2 (t) + Y 2 (t)*, 

with Y 2 (t) being an N x N matrix with all entries independent complex standard Brow- 
^ nian motions. 

X 2 (t) takes values in the complex hermitian matrices, and is distributed like a scaled 

GUE matrix for fixed time t. Using the Ito formula and the invariance of the increments 

under unitary transformations, it is possible to prove that the law in path space of 
CNJ the eigenvalues X\{t), Ajv(t) of X 2 (t) is the unique weak solution of the system of 

stochastic differential equations (SDE) 

8 P) dA lW =A dwVt) + 2 ^_i_ d , fc=l JV 

£f-} l^k 

with Wk independent standard Brownian motions and j3 = 2 in our case. Although it 
is not hard to formally derive ([2]), the singular level repulsion terms (which originate 
from second order perturbation theory) cause that a solution of (|2j is only defined up 
to the first collision of two eigenvalues. Thus, to make the derivation of (|2j from 
rigorous, one needs to prove that almost surely no level crossing can occur. This was 
^ accomplished directly by analysis of |TJ in [17J, and in |7 ) ,[[2~T | via a Lyapunov function 

argument for (|2j (for general f3 > 1), which is presented comprehensively in |2|. 

Among other applications in physics, random matrices have been employed to un- 
derstand the properties of disordered conductors, in particular the phenomenon of uni- 
versal conductance fluctuations ID EH- In the arguably simplest case, one can model 
the conductor as a quasi-lD wire, a system of iV interacting channels. The matrix in 
consideration is the transfer matrix M g ^Nx2N w ]- l j c ]- l ma p S the quantum state at the 
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left side of the wire to the state at the right side (rather than mapping incoming states 
to outgoing states, as a scattering matrix does). 

Instead of the self-adjointness of Wigner matrices, the main constraint for transfer 
matrices is current conservation, (corresponding to unitarity of the scattering matrix), 
which in a suitable basis reads 

(3) M*Z Z M = E z , S z = 

ffl9| , and we define the group 

(4) g 2 = {Me C 2Nx2N : M*Z Z M = S 2 } . 

In case the underlying quantum mechanical system is time-reversal symmetric, i.e. in the 
absence of magnetic fields, the transfer matrix additionally has to satisfy 

(5) S.jWE, =M, E x = 
Accordingly, 

(6) Gi = {M£ C 2Nx2N : M*Y, Z M = S 2 and = M} . 

In addition to the groups Qr with parameter /3 = 1, 2 as defined above, the value (3 = 4 
is also algebraically and physically meaningful, as it corresponds to a quantum system 
allowing for spin-orbit scattering, (5|. However, this last case does not introduce any 
new stochastic features and is disregarded for the sake of notational simplicity. One 
should remark that the correspondence of symmetries of the quantum system and the 
choice of matrix ensembles already arises for Gaussian ensembles and Dyson Brownian 
motion lH"8l , the ensemble suitable for time-reversal invariant quantum systems (and 
thus losely related to our Q\) being 

(7) X 1 (t) = Y l (t) + Y l (tf, 

with Y\ (t) being an N x N matrix with all entries independent real standard Brownian 
motions. 

As in the theory of Wigner matrices, the focus is on the spectral properties of those 
transfer matrices: for each M. there exists a decomposition 

M _(M ++ M + \_(U + 0\( A (A 2 -l)^A^ 0\ 

(8) M -{m- + M-)-[o u„){ {a ._^ a ){o V-) 

with (non-unique) N x N unitary matrices U + , U- or V + , V_, and a diagonal matrix 
A > In- For general M. G Q2, there is no relationship between those unitaries, while for 
M G Q\, U + = U- and V+ = V-. For later use we define the submanifolds Qp as the set 
of those matrices in Qr with a decomposition such that A > ljy and has nondegenerate 
eigenvalues. 

For physically realistic conductors, M. will be a random quantity due to microscopic 
disorder, so that on a macroscopic level it is natural to model the transfer matrix as a 
stochastic process {M(s)) s>0 depending on the length s of the wire. 
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If we combine two pieces of wire with transfer matrices M.\ and M.\i, respectively, 
we expect the concatenated wire to have 

(9) M = MuMi 

as its transfer matrix [19J and if we furthermore assume that the transfer matrices of 
disjoint short pieces of wire are i.i.d., we arrive at the stochastic differential equation 

(10) dM(s) = dC(s)M(s) (s > 0). 

with a suitable stochastic process (jC(s)) s>0 that encodes the scattering properties of 
short wire pieces and is chosen such that (|3f (and, for f3 = 1, ^) are conserved. Stochas- 
tically, every intial distribution for pi)} makes sense, but it is physically reasonable to 
assign to a zero-length conductor the transfer matrix 

(11) M(0) = h N - 

The process C we will actually consider is a 2N x 2N matrix-valued Brownian motion 

«+W &M 



(12) • *•(.) M » 

with N x N block processes o + , a_, b which, in terms of ( |1|7| , are distributed like 

a+ = a_ =iX 2 /V2N 

(13) Reb = Imb = X 1 /(2 v / FTl) (0 = 1) 

b = y 2 /Vn ((3 = 2). 

For /3 = 1, Reb and Imb are mutually independent and jointly independent of a±, while 
a + and o_ are correlated by a_ = o+. For (3 = 2, however, the processes a + , a_ and b 
are all independent. 

One motivation for this choice of C is that it predicts the same dynamics for the di- 
agonal matrix A as the maximum entropy assumption MT9| . The reason is, that after 
endowing the groups Qp with a suitable right-invariant metric, ( [To] ) is just the Brown- 
ian motion on the Riemannian manifold Qp, while ||T9l essentially single out the heat 



kernel by their maximum entropy requirement for the propagator. Moreover, (10 1 has 
been derived in 0111 as scaling limit from a microscopic quantum model (an Anderson 
model on a tube) — up to minor deviations originating from the Anderson Hamiltonian 
dictating a preferred basis in the weak-perturbation limit. 

Finally, and most important in the present context, the increments of C are invariant 
under the conjugation 

(14) WdCU = dC. 

U can be any unitary matrix of the form 

'U+ 



(15) U- . t Q v _ 
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with the blocks U+ and U- chosen independently for f3 = 2, and correlated by U + = U- 
in case (3 = 1, compare (21), (22) in [3J. This perfectly corresponds to the orthogo- 
nal/ unitary invariance of the GUE/ GOE-like increments of ( l]7!, and will finally allow 
us to prove an autonomous (i.e. eigenvector-independent) equation for A. 

First, we give a formal derivation and start with a stochastic evolution for the ma- 
trix M* ++ M++ and its eigenvalues X k , ordered like Ai > ... > Xjy. Inserting the ++ 
component of pO) , 

(16) dM++ = daM++ + dbM-+, 
into Ito's formula, we obtain 

d(M* ++ M++) =M* ++ (da+M++ + dbM-+) 

+ (M* ++ da* + + Ml+db*) M++ 

(17) w 

+ (M* ++ da* + + M*_ + db*) (da+M++ + dbM-+) 

=M\ + dbM-+ + M*_ + db*M ++ + M* ++ M ++ ds + M*_ + M-+ds 

where we used the independence of o + and b, and their explicit form ( [13] ) in the last 
line. By second order perturbation theory for the eigenvalues of M* ++ M++ and Ito's 
formula, 

(18) d\ k = (v k ,d(M* ++ M++)v k )+J2 

l^k 

with v k the eigenvectors of M* ++ M ++ . In view of |8|, v k are the columns of V*, while 
-y/Afc = A kk are the singular values of M.++. Furthermore, by 

(19) uim- + v; = {k 2 -iy 

U* + M ++ V* = A 
it is convenient to consider the process with 

(20) db(s) = Ul(s)db{s)U-(s) 
for which, by p4|, 



[vi,d(M* ++ M ++ )v k ) 
Afc — Xi 



(21) db = db. 
With these definitions, 

dA fc =2 v / X k {X k - l)d (Reb fcfc ) + (2A fc - l)ds 

(22) + (VMAfc _ l)dbjjfc + A/A fc (Ai - l)db M 



(V^(A fe -l)db/ fe + VA fc (Aj - l)db M ) /(A* - AO- 



We can use (21 1 with (T13l to see 



(23) Reb fcfc (s) = -J p {N \ ) + 2 Bk(s) 
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for all k = 1, N, with independent standard real Brownian motions B k , while 

(VAKA fe - l)db ifc + v%(Az - Tjdbiaj (y/Wk - l)db^ + v%(Af - l)db w 
( 24 ) /3 ^ 



/3(iV-l) + 2 

and thus 



(2A fc A; - A fc - A;)ds, 



r>*\ *\ /ox P 2A fc A ' ~ Afc - A; \ 4A fc (A fc -l) ^ p 

(25) = " 1 + W-D + 2 g A, -A, J ds " V m -l) + 2 dB - 

After introducing the transmission eigenvalues T k = A^ 1 G [0, 1], another application of 
the Ito formula yields the DMPK (Dorokhov-Mello-Pereyra-Kumar) equation for the 
transmission eigenvalues (compare (3.9) in [5[, (24) in Il3l): 

(26) dT k (s) = v k (T(s))ds + D k (T(s))dB k (s), 
B k , k = 1, N independent Brownian motions, 



2T k / 8 T k + Tj - 2T k Tj 



(27) 



/3iV + 2 - /3 2 ^ T k — Tj 



Tl(\-T k 



y /3iV + 2-/3' 

Remark. As already noted in [3], the term "DMPK equation" usually refers to the for- 
ward equation for the density of the T k 's. We also use it for the corresponding stochastic 
differential equation, as this is the more natural object in our analysis. 
The form of both the drift term v k and the diffusion term D k causes the transmission 
values to decay to as s tends to infinity, which perfectly matches the decrease of con- 
ductance with increasing wire length. However, the drift v k also contains repulsion 
terms originating from second order perturbation theory. As a consequence, the eigen- 
values T k "try to avoid" degeneracy. Thus, and similar to Dyson Brownian motion, 
the above derivation is merely "formal", as the Ito formula is only applicable if the de- 
nominator \ k — Xi (or T k — T{) stays away from zero, i.e. M* ++ (s)M.J t +(s) never has 
degenerate eigenvalues. For s > 0, this is already a nontrivial finding, but even more, 
we want to start the evolution equation pD) with the completely degenerate matrix pTj l. 
The goal of our paper is the derivation of the following rigorous statement about the 
transmission value process of (M(s)) s>0 . 



Theorem 1. Let 8 = 1 or B = 2, and (M(s)) s>0 be the solution of (10), starting from (11 1. 
Then the distribution in pathspace of its transmission eigenvalue process is given by the law of 
the unique continuous process {T k (s)) s>0 k=1 N which starts from T k (0) = lfor k = 1, N 
and is a strong solution to ( (26| for s > 0. Moreover, the transmission eigenvalues T k , Ti for 
k / I almost surely never collide for s > 0, with each staying in the open interval (0, 1). 

Speaking in terms of manifolds, we thus see that the paths of the Brownian motion 



(101 on Qp, even though starting from a degenerate matrix M(0) € Gp\ Gp, are almost 



surely contained in Gp for all positive s, and do not explode for finite s (this would 



6 



M. BUTZ 



correspond to T k (s) = 0). It will be obvious from the proof of Theorem |4]that the 
way to handle the singular initial condition ( TIT) could as well be applied to any other 
degenerate M(0) G Gp\Gp- For the sake of notational simplicity we stay with the 
physically most interesting case M(0) = l2N- 

We will start from considering the DMPK equation not as an eigenvalue process, but 
a process in its own right and show that it has a unique strong solution for all /3 > 1 in 



subsection 2.1 Once this is established, we verify in subsection 2.2 that the transmission 
eigenvalue process of M has the same law in path space. 

Whenever possible, the analysis will be carried out in the picture of the transmission 
eigenvalues T k (rather than the \ k ). This is because the T k are the physically most 
intuitive quantity, each T k representing the transmittance of one channel. Accordingly, 
when modelling the passage of electrons through a iV-channel wire, the dimensionless 
conductance of the wire is given by the Landauer formula (|5|, equation (33)), 

N 

(28) g = Y. T *- 

k=l 

This formula for g and other linear functionals of the T k describe the most important 
transport properties of the disordered conductor. Plugging the DMPK forward equa- 
tion into ( |28| l, [6j established the universal conductance fluctuations 

(29) Var( 5 ) = ~ 

for P = 2 in the intermediate metallic phase with 1 <C s <C N. 

2. Rigorous analysis of the process of transmission eigenvalues 

2.1. The unique strong solution of the DMPK equation. We first investigate equation 
( p6| starting from a nondegenerate initial condition < T\(0) < ... < T/v(0) < 1, 
corresponding to .M(O) G Gp- 

Theorem 2. Let (O, P) be a probability space endowed with a filtration T = (Ft) t>0 and 
Bi, Bn be independent real Brownian motions adapted to T. Let Dn be the set of strictly 
ordered N-tupels in (0, 1), = {x G M. : < x\ < ... < xn < 1}- Assume initial condi- 
tions T(0) = (Ti(0),...,Tjv(0)) G D N , and take f3 > I. Then the stochastic differential 



equation (26* has a unique strong solution with initial condition T(0), and T(s) G Dn for all 



s > 0. The law of this strong solution also defines the unique weak solution o/([26|. 

Proof. Following the ideas of [2]. For R > 1, we regularize the singular repulsion in the 
drift term of |27| | by setting 

(30) Vk {T) = XR{T) [- Tk + WT2-^[ l - Tk+ 2^ k T k -T 3 



d[ R \t) = XR (T)D k (T) 



with xr a smooth cut-off function compactly supported in the interior of Dn and 

XR (T) = 1 whenever T G D N and dist(T,dD N ) > (V2R)-\ Now v[ R) and d[ R) 
are both uniformly Lipschitz, so the modified system has a unique strong solution 
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(T( R ) ( s )) s>0 an d a unique weak solution (|15J, Theorems 5.2.5, 5.2.9, |l2"0l| . Lemma 5.3.1). 
he idea of the proof is that T^ RS) (s) should be a solution of the original DMPK equation 



(26 1 as long as the differences of the 1): and their distance to the boundaries of [0, 1] 
are larger than R . Taking R to zero should then yield a solution to (26 for all times. 



To control the separation of the eigenvalues and their distance to and 1, we define a 
Lyapunov function 



N I N 

(3i) f(T) = J2 -2iog(|r fc |)-2io g (|i-r fc |)- Yl Ml?* -ill 

k=l \ l=X,l^k 

It will be enough to define / for T from D^, as we will always use it together with 
stopping times which will make sure T doesn't leave this set. 



Note that if we take T € D^, all the logarithms in (31 1 are positive, thus we have 
f(T) > and also 

(32) -log(|T fc -T,|)</(r). 
Define 

(33) S M ,R = inf {s > : / (t^ (a)) > m} 

for M > 0. As / is smooth on sets where it is uniformly bounded, Sm,r is a stopping 
time. On the event {Sm,r > s}, we have 

(34) - log(|Tf \s) - T{ R \s)\) < f(T^(s)) < M 
so if we take M = M(R) = logR, 

\Tt R) (s)-Tl R \s)\>e- M = R-i 

(35) \T [ k R \s)\>e- M = R- 1 

\l-T {R \s)\ >e~ M = R-\ 

Thus, up to S M ( R ) R , T( R) and all with R > R coincide and T^ R) also solves the 
original equation d26| up to Sm(r),r- 

Therefore, to calculate df(s) up to the stopping time Sm = Sm(r),r it is enough to 



plug (26 ' into Ito's formula (ffl5|. Theorem 3.3.3). We drop the overscript (R), as it does 
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not matter in this context. 



(36) 

d/(T( S )) = 

22J ^ + i is* r *- 2 i- 



fc=i 



f 2T fc / A T k + Ti- 2T k Ti . , 



+ s(i + ( i - 1 ^ 2+ ^ (rfc -^ )2 ) 

+ d,4(s), 

with (yl(s)) s>0 denoting the local martingale fulfilling 



k=l \ j=lj^k 



(Remark: Because of ([35]), A(s A Sm) is certainly a martingale, which is enough for 
our purpose. We only see that A is really a local martingale after having verified that 
Sm /* oo as M — > oo.) 

Introducing Tq = 0, T^+i = 1, we can more conveniently rewrite df 



N / N+l 



BN + 2 - B ^ \ ^ T k -Ti 

H k=l \j=0,j^k K 



N+l 

, (B - 1)T| + BNT k (T k -l) + B V - Pk - I ds 
(38) \ {P ' k ' ,4^-7] I 



x 



+ d4(«), 
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where we have set pt = T%(1 — Using po = Pjv+i = 0, we have 

N ( ( 1 \ ( Pk \ Pk 

Z ~ E^ TV _ T, E^ T, - T, ~~ E^ 



fc=l \ \j=0,j^k Tk Tj J \l=0,l^k Tk Tl J m =0,m^k Tm " > 
N+l / / N+l 1 \ / AT+1 ^ \ AT+1 

EZ I I E^ j 1 , — T ) I E^ Ti — Ti I ~~ ^ fTi, — T ) 2 

fc=0 \ \j=0J^fc K J / \l=0,l^k ft 7 m =0,m^fc V ' 
N+l N+1 N+1 

Pk 



(39) 



T, Tk — TiJ Tj — Ti 



N+1 N+1 N+1 , 

fc=0 j=0 z=o ^ k 

N+1 N+1 N+1 , x 

fe=0 i=0 Z=0 \ K j K 1/ j i 

N+1N+1N+1 

fc=o i=o i=o 
and thus for T G L>jv 

(40) |Z|<?(iV + l) 3 . 

Furthermore 
(41) 

f ^ T fc( l-T fc) ^ Ml-Tj-Tfr-iy 1 V 

2. 2. Tfc _ r -l^l^ 1 ^^^ 2 (T fc - t 7 ) aZJ 1 "^- 1 ^ 

k=lj=0,j^k K 3 fc=0 fc=0 v K J ' fc^j 



SO 

N N+l 

(42) 



E E 

k=l j=0,j^k 



T fc (l-T fc ) 



(N + l)(N + 2) 



for all T G D N . Finally, 



N N+l „ 2 Af W _ 2 

E V __k__ _ _ , _z_k__ , \ ^ J fc 

^ T,- - T fc ^ fe ^ 1 - T fc ^ T,- — T fc 

(43) k=ij=o,j^k 3 fe fc=i fc=i fc j,k=i,&k 3 

^ ' N N Ar 



\ \ ^ T fc + T i , \ ^ T fc 

2 ^ 1 - T* 

fe=i i=i fc=i * 

and thus on D 

JV N+l T 2 N T 2 

(44, -^< E E ^< E _5_. 

fe=l j=0,j^k 3 K k=l K 
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Thus there is a continous bounded function g with g(T) < CN 2 for T G such 
that 



d/(T( S )) 



4(0-1) " / T| 



(45) =^->» + i^ 




Recalling that /3 > 1 and A (s A 5m ) is a martingale, 

(46) E(/(T( S A S M ))) < CN 2 s + /(T(0)). 
As / is nonnegative, we have the Markov estimate 

(47) P (5 M < «) < jjW(T(s A 5 M ))) < ° N2s + J im) 

for all s > and M > 0. 
Now as by definition 

(48) 5 M = S M ^x > 5 M ( jR ) )iJ = 5m 

for R> R,we conclude from P (S M < s) -> as M oo for all s > that 

(49) 5 M / oo (M -> oo) 



almost surely. As we have already seen that T( R > is a solution to (26 1 up to S M ( R ) Rr 
setting T = T^ R ' up to Sm(r),r for all R > 1 provides us with a strong solution of (26 
for s G [0, oo). This finishes the strong existence part. 

To prove strong uniqueness, let a strong solution T of (26 1 be given and define the 
stopping time 5m quite analogously to 5m, but this time based on T instead of T^ R \ 

(50) S M = inf {s > : / (T(s)) > M} . 



If we take e M = R, we have a result like (35 >, and by strong uniqueness of T^ R \ T and 
T( R ) coincide up to 5m- As one can show 

(51) 5m / oo (M — >• oo) 

just as before, T has to be the solution we obtained in the existence part almost surely 
for all times. 

Weak uniqueness of T follows by exactly the same reasoning from the weak unique- 
ness property of for all R > 1. 

Finally, T(s) £ for all s > almost surely, because we have seen that for any 
s > there is almost surely a (random) R > 1 such that T(i) = T^ R \t) for all t < s, 
but the stay in D^r almost surely due to the cutoff in the definitions of and 
Df\ ' □ 
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Next, we want to extend our existence and uniqueness result to the initial condition 
that corresponds to M(0) = I2N, i-e. T k (0) = 1 for all k = 1, N. To do so, we adapt 
Proposition 4.3.5 and Lemma 4.3.6 about Dyson Brownian Motion in [2] to the DMPK 
equation p6| . 

To control the convergence of approximate solutions, we start with a comparison 
result for solutions with different initial conditions in D^. 



Lemma 3. Fix (3 > 1, N E N and let (R(s)) s>0 and {T(s)) s>0 be two strong solutions to (26* 
with i?(0),T(0) g D N and 

(52) Rj (0) > Tj (0) for all j = 1 , . . . , N. 
Then we almost surely have for all s > and all j = 1,...,N 

(53) Rj(s) > Tj(s) for all j = 1, ...,N. 
Proof. We start with defining the stopping time 

(54) S e = inf js > : minmin{l - Rj(s), 1 - Tj(s)} < e 

for e > 0. As we have seen in the proof of Theorem |2j 

(55) S e — > 00 a.s. 
as e — > 0. 

Furthermore let the stopping time S be when the ordering |53] l is violated first, i.e. 

(56) S = inf {s > : Rj (s) = Tj (s) for some j} . 

We want to show that S = 00 a.s. 

From (27 1 we see that Y^k v k(T) is a polynomial in the Tj/s so that 



(57) 



5> fc (i?( S ))-5> fc (T( a )) 



A' 



iY 



<C7^|« fc ( S ) 

fc=i fc=i 

on the event {S > s}, where we have used that all Tj, Rj are contained in the bounded 
set [0, 1]. Here and in subsequent estimates we use constants C, C... which might de- 
pend on N, but not on e. On the other hand, for 



(58) 



D k (T k 



T 2 Al~T k ) 



/3N + 2-/3 

we have \D k (R k ) - D k (T k )\ <C\R k - T k \ jjl as long as T k , R k < 1 - e, so 

N 



(59) 



( D k (R k (s))-D k (T k (s))Y 



< ^=Y^(R k (s)-T k (s 
v e k=i 



)) 



on {S A S € > s}. 
Thus 

(60) 



iV 



u(s) = Y J (Rk(s)-T k (s)) 



k=i 
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solves the stochastic differential equation 

(61) du(s) = f{R{s),T(s))ds + a(R(s),T(s))dB(s) 
with B a Brownian motion and 

\f(R(s),T(s))\<Cu 

(62) C 

\a(R(s),T( S ))\<-^u 

on {S A S e > s}. By Ito's rule, 

1 



f a(R,T)dB > -oo 



(63) d Iog(«(a)) = -^-du(s) ^a 2 ds = h(R, T)ds + a(R, T)dB(s) 

u{s) 2u{s)' L 

with h > -C - C 2 /e and \a\ < C/y/e. Thus, on the event {S A S € > s}, 

(64) log(«(a)) > log(n(0)) - (C + C 2 /e)s + f er(fl(t), T(i))dB(i). 

Jo 

As 

rSAS e As fi2 B 

(65) / a 2 dt < — , 

we have 
(66) 

almost surely on {5 A 5 e > s}, and thus from ( |52| 

(67) u{s) > u(0) exp f-(C + C 2 /e)s + a(i?, T)dB^j > 

almost surely on {S A S € > s} for all e > 0. So with S e — > oo as e — ► 0, u(s) > almost 
surely whenever s < S. 

Now assume that for some finite so, IP (S < so) > 0. And define Uj(s) = Rj(s)—Tj(s), 
which is almost surely nonnegative for all s < S by continuity of paths. Whenever 
S < oo, there is a (random) fc G {1, N} such that Uk(S) = 0, but as we have seen 
above that u(s) > almost surely on {s < S}, there is almost surely a j such that 
R j (S)-T j ( i S)>0. 

For all k, Uk solves the stochastic differential equation 

(68) 

du k (s) 

>l A " ' ] M s ) ~ o,* T 2 ^ , J Rk(.s) + T k (s))u k (s) + r(R(s),T(s)) ) ds 



/9(JV — 1) + 2 v ' /3(JV — 1) + 2 



, Rl(l-R h ) T 2 (l-T k ) T> 
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with r denoting the difference of the repulsion terms 
(69) 

r(R,T) 



'pN + 2- 



/? R k + Rj-2R k Rj \ _ Tk fT k + T 5 - 2T k T :j 



Rk — Rj J V T k — Tj 



(3 ^--^ RkT k (Rk — T k ) + 2TjRj(R\ — T|) + (R k + T k )T k (Rj — Tj 



(3N + 2-(3^ (R k -Rj)(T k -T, 



3i 



P \ Tj (T k + R k + Rj + 2R k T k ) (T k — R k ) 

+ l5N + 2 -pZ^ (R k -Rj)(T k -Tj) 

s dependence suppressed. In the last two lines, the denominator (R k — Rj)(T k — Tj) 
is always positive as R, T G D^. Now on a path with S < oo, s < S, and k such 
that u k (S) = 0, the numerator in the second last line is positive, the one in the last 
line is negative. However, while the negative factor (T k — R k )(s) = —u k (s) vanishes as 
s /* S, some of the positive terms Rj — Tj are bounded away from zero uniformly in a 
neighborhood of 5, so there is an (w-dependent) s\ < S such that 

(70) r(R(s),T(s)) > for all s G [s±, S] . 

Then given our assumption that P (S(u) < so) > 0, we can find deterministic, constant 
s\,e and k, such that the event 

(71) A = {s x < S < s } n {r(R(s),T(s)) > OVs G [s u S]} D {u k (S) = 0} n {S e > s } 
has a positive probability P (A) > 0. On the event A we have for s G [s\, S] 

(72) du fe (s) > -2u k (s)ds + a k (R k (s),T k (s)) dB k (s) 
with 

(73) \a k (R k (s),T k (s))\<-^u k (s). 
As before we obtain 

(74) u fc (s) > u fc (ai) exp ^-(2 + C 2 /e)(s - si) + ^ *fcdB fc ^ 
for all s G [si, 5]. Again, we have <t| < C 2 /e, and thus 

(75) / a k dB k > -oo 

for almost surely on A. Together with u k (s\) > on A this means for u k (S) > on A, 
contradicting the assumption u k (S) = 0. So our initial assumption is wrong and S(co) 
has to be infinite for almost all u; G O. This means that the initial ordering |52] l stays 
conserved almost surely. □ 

Theorem 4. Let f3 > 1, N G N, and Tjt(0) = 1, = l,...,iV. Tfoen fere zs a unique 
continuous (T(s)) s>0 such that T(s) G Dn for all s > and (T(s)) s>0 zs a strong solution 
to the DMPK equation ( |26] >. Again, the law of this strong solution also defines the unique weak 
solution. 
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Proof. To approximate T(0) by intial conditions T^ n >(0) £% define 



N + l-k 



n 



those initial conditions allow for a solution (T^ (s)) s>0 
\n, s. By Lemma pi 



Tt\s)<T^ (s) 



(76) 

for all n > N+l. By Theorem 
of g with rW(s) G Dn for a 

(77) 

almost surely for all N + 1 < m < n, s > 0, so the converge from below to a limit 
process (T(s)) s>0 with T fc (0) = 1 for all k and T k (s) G for all s > 0. 

This convergence result can be considerably strengthened. Define for m < n 

N N 

(78) u^ n \t) = \T [ k\t) ~ T, (m) (t)| = ]T (Tl n \t) - T^ m) (t)) > 0, 



k=l 



k=l 



which solves the stochastic differential equation 

(79) du( m - n ) (t) = / M m ) (t) , (t)) dt + a M m > (t) , (t)\ ab^ 1 ^ (t) 

with £?( m > n ) a Brownian motion. This time we estimate 

/ (T^ (t) , (i)) < Cu^ (t) 

T^ m \t),T^ n \t) 



(80) 



for all i, and only iV-dependent constants C, C. Then the stochastic process 
(81) 

,(m,n) (q\ _ 



is a continuous martingale and its square Y m ^ n = n is a continuous submartingale. 

As u < AT, CT m '"(t) = a {T^ m \t),T^ n \t)) is bounded, and as T^ m) (s) -> T fc (s) from 
below for all fc, s as m — > oo, we also have 



(82) 



sup |a m ' n (t)| ^0 (m^oo) 



n>m 



almost surely and for all t > 0. Using this and dominated convergence in the Ito isom- 
etry we obtain 

(83) K m = sup E (Y m , n (s)) = sup / E (\a m > n (t)\ 2 ) dt 

n>m n>mJo ^ ' 

as m — ^ OO. 

The key idea for uniform convergence is Theorem 1.3.8 in ||T5l , which states that for 

any a > 0, s > 

/o/i\ to I Ay ("jA \ I ^ (*m,n( s )) K m 

(84) sup P I sup r m ,n(*J > Q! I < sup ! = — . 

n>m \0<t<s ' J n>m OL a 
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Choosing ej > and a subsequence rrij — > oo such that 

(85) f j (^L + V ^) <(X 

we obtain 

(86) > 'P ( sup Y m ^ mi+1 (t) >£, ]<> ^ < oo. 



i=l ^ ej 



CO / \ CO 

.7 = 1 .7 = 1 J 



By the Borel-Cantelli lemma we thus have an almost surely finite J(co) such that for all 

3>J 

(87) sup Y m . mj+l (t) < ej 

0<t<s 

and thus for all t < s, j > J 

u N.^+i)(f) < /^- + n (^Wi)(o) + C / «^' m ^+ 1 )(r)dT 

(88) ' -A) 

u^'^'fi) < (^7 + u (m ^' m ^+ l) (0)) e CT , 

where we used the Gronwall inequality. With the monotonicity of the in m and 
thus of almost surely have 



lim sup sup u {m ' n \t)= lim sup r«( m i' m «'(t) 

m ^°°n>mo<t<s o<t<s rrf 

— — J>J 



(89) 



N 

< lim e^(E^+E ( T ^°) " T t l) (°) 







and thus almost surely uniform convergence of 2^ to on compact intervals for 
all k. (This means the speed of convergence is uniform in compact intervals [0, s], but 
might still depend on the particular path). In particular, all are continuous almost 
surely. 

Now assume for a moment that we already know that for a certain so > 0, T(sq) € 



Dn almost surely. Then denote by ( T(t) ) the unique strong solution of (26 starting 

V / t>so 

at so with T(so) = T(sq), which exists by Theorem[2j 
Define S e = inf js > s : f N (s) > 1 - ej for e > and 

N 

(90) u^(t) = J2(fk(t)-Tj i m) (t)) 

we have functions /, a and a Brownian motion i?( m ) such that 

(91) du (m) (s) = / (f(s),T( m \s)) ds + a (f(s),T^(s)) d£ (m) (s) 
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with the same estimates as in d62l 



(92) 



/(f( s ),r( 



a 



on the event |so < s < S e \. From these estimates and E (u(' m )(so) 2 ) 
one can conclude by the standard Gronwall method that 



as m — > oo, 



(93) 



E f (sAS e y 







as m — > for all s > sq and all e > 0. By the monotonicity of in m we therefore also 
have u( m ^ (^s A 5 £ J — )■ almost surely for all s > sq and e > 0, and as S e /• almost 

surely as e — >■ 0, even u( m ) (s) — >■ almost surely for all s. Thus (T(t)) t> and ( T(t) ) 

are modifications of each other, and by their continuity they are indistinguishable. 

To verify that (T(s)) s>0 is a solution to (261, the only thing left to prove is that there 
is almost surely no si > such that T(s) £ Djy \ Dn for s G [0, si). First, almost 

surely T fc (s) > r^ 1} > for all s > 0, k = 1, iV, by Theorem^and Lemma|3| Next, 
assume < Tk(s) < 1 for all s > 0, and all k, but that with positive probability there is 
a neighborhood of zero where some of the components of T "touch". Then, due to the 
continuity of T as a locally uniform limit of continuous functions, it is possible to find 
a deterministic k, a fixed open interval (a, b) C M + and a fixed 5 > such that on an 
event ft' with P (ft') > for all s G (a,b),T k (s) = T k+1 (s), but T k+1 (s) < T k+l (s) - 5 for 
all I > 2, and Tn(s) < 1 — d. On this event ft', with the uniform convergence of T'"' to 
T, we have for sufficiently large n, for all s G (a, 6) 



(94) 



(n) 



k+1 



/(rM)d s + a(TMW +1 (s 



/ (r (n) ) 



> Ci 



c-l 



1) 



with a uniformly bounded but / becoming more and more singular because of the 
and repulsion 

(95) 

with 5 fixed while e n 
(96) 

which is not possible for a solution of p6| with intial conditions in Dn due to Theorem 

m ' 

Finally we have to exclude that there is a neighborhood of such that Tn(s) = 1 for 
all s G [0, si). If so, there was an open interval (a, b) again, a k and a 5 > such that 
on a positive probability set ft", T k (s) = 1 for all s G (a, 6), but Tj < 1 — 5 for all j < k. 



0. Thus, with probability P (ft') > 0, there is an n G N such that 

rW(6) > i 
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Then almost surely in Q," we have the following convergences in 
(97) 

-,(n) / a rp{ri) , T {n) _ 0T (n) T (n)' 



T( n A (A < -T (n) + k 1 - T (n) 4 " V" k j k j 

: V ) [S) - ±k ^f5N + 2-(3\ k T (n)_ T (n) 

p{k-\) 



< 



f3N + 2 - (3 
2 







pN+2-p 

and 

(98) D fc (t^) (s) 
for all s G (a, 6) as n — >■ oo. By Ito's isometry and dominated convergence, 

(99) f D k (TW) (s)dB k (s) 

•/ ft 

in L 2 (S7", P) as n — >■ and thus almost surely on a subsequence rij, but then 

(100) H(.) = .lim !**>(.) < 1 - ^^^(» - «) 

almost surely on f2" for all s G (a, 6), contradicting our assumption. 

Altogether, (T(s)) s>0 is almost surely continuous with the right initial conditions, and 

(T(s)) s>Q is a strong solution of (|26]l. 



To finish the proof of Theorem |4j we have to see why T is the unique strong solution 
to (26 i with the given initial conditions. Assume (T(s)) s>Q was another continuous 
process solving (26 1 for s > and starting with T k (0) = 1 for all k. For fixed n, as 

1^(0) < 1 for all k, and by continuity of T and T^™), there is a random, but almost 
surely positive S such that still 

(101) T( n \s)<T k (s) 

for fe = 1, N on the event {S > s}. As T and are solutions of (26 i for all positive 

s, we apply Lemma [3] to obtain T^ n \s) < T k (s) almost surely for all s > 0, and taking 
the limit in n, also 

(102) T k (s) < T fc (s) 

for all and s. This means we can use the random variable 

N N 

(103) u = Y J \Tk-T k \=Y,{T k -T k ) 

k=l k=l 

to control the distance of the two solutions. We have 



(104) 



du(t) = f {T(t),T(t)) dt + a [T(t),T(t)) dB(t) 
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for all t > 0, with B a Brownian motion and 

\f{T(t),T(t))\<Cu(t) 

(105) _ , _ , 

\a{T(t),T(t))\<cVW) 

with iV-dependent constants C, C. We pick an s > and estimate 

(106) E(u(t))= rE(/(T(r),r(r)))dr<E(n( S ))exp(C7(t- S )) 

is 

for all t > s. Now as u < iV, and u(s) — )• as s — > due to continuity of T and T, we 
have by dominated convergence in ( 106| 

(107) E(|u(i)|) = 

for all t > 0, and again by continuity of T and T, T and T are undistinguishable. Thus 
we have obtained the unique strong solution by our construction. 

For the proof of weak uniqueness, we note that there is an equivalent of Lemma [3] 
that allows us to compare a weak and a strong solution with initial data ordered like 
(52 1, which allows us to repeat the above argument with (T, B) a weak solution of the 
DMPK equation for positive times, while T^(0) = 1 for all k. □ 

2.2. Equality in distribution. Now we are ready to prove that the DMPK equation 
actually describes the process of the transmission eigenvalues. 

Let /3 = 1 or j3 = 2, and (0,P, J 7 ) be a probability space with a filtration {J : s ) s >q, 
and, using Brownian motions with respect to this filtration define N x N matrix valued 
processes o + , a_, b and o + , a_, b distributed according to (13 ', with (o+, cu, b) indepen- 
dent of (o+, a_, b), and both triples obeying the dependence /independence structure 
described under equation ( [T3| . 

As it is convenient for later calculations, note the component-wise distribution of the 
matrix blocks: a + , a_, d + and d_ are all equal in distribution with \w entries 

t l/x/2JV(^( S ) + ^( S )) v<v 

(108) ^ u (s) = l i/VNjB^js)) ^ = ^ 

-a Ufl (s) n> v 

for fj,, v = 1, N, with all components ofB R and B 1 independent, standard real Brow- 
nian motions. The blocks b and b are equal in distribution with entries distributed like 

t l/y/2(N + l)(B*(s) + i% v {s)) fl<U 

(109) b, u {s) = y/l/(N + l)(B^(s) + iB^is)) fi = fi (/3 = 1) 

I b Ufl (s) fj,> v 

or 

(110) b^s) = l/V2N{B*( a ) + iB^(s)) (P = 2), 

again with all components of B R and B 1 independent, standard real Brownian motions. 
Define (C(s)) s > and (M(s)) s > by 

r( v ( a+{s) b(s) 
£(S)= U*( S ) o_W 
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and 

(112) dM(s) = dC(s)M(s) 

with initial condition M(0) = hN- Let X M (s) G 1^ be eigenvalues of M* ++ M ++ {s), 
ordered by A^ > ... > A^ 1 and define the transmission eigenvalues as before by T^ 4 = 

(Afr 1 . 

On the other hand, define N independent Brownian motions by p3| , so 

(113) B k (s) = - v / /3(iV-l) + 2Reb fcfc ( S ), 



and let (T k (s)) s>0 , k = 1, ...,N be the solution of the DMPK equation (26 driven by 
those Bk, starting from Tfc(O) = 1 for all k, as constructed in TheoremE] 



Theorem 5. The distributions of (T M (s)) s>Q and (T(s)) s>0 in (C r n[0, oo),B (C r jv[0, oo))) 
are identical. 

Remark. The Borel c-field B {C r n [0, oo)) is defined with respect to the metric of locally 
uniform convergence. For the statement of the theorem, we obviously do not need T M 
and T to be defined on the same probability space, but this will be very convenient for 
the proof. 

The idea of the proof is similar to that for the Dyson Brownian Motion in [2J, es- 
sentially by reading the formal derivation of ( 26 1 backwards and making every step 
rigorous. To get started, note that we know that T M (s) G D N , but even have 

Lemma 6. For T M defined as above, 

(114) limP(r- M (s) G D N ) = 1. 



> , 



Remark. Once we have shown Theorem |5j we will actually know much more, namely 

(115) P({T M (s) G DjvVs > 0}) = 1. 

For the moment, we morally need P (T M (s) G £>at) = 1 for any fixed time s > 0. This 
could be achieved by writing ( 112| l as a Brownian motion on the manifold of transfer 
matrices, and then using the smoothness of the heat kernel (3 [131 to verify that the 
lower-dimensional manifold corresponding to T M (s) G Dn \ Djq has probability zero. 
Instead, we use the weaker statement ( 114| which can be proven directly and suffices 
as well. 

Proof. T M (s) G Djv is equivalent to \ M G (1,00)^ and nondegenerate, and finally, 

that the matrix A4*L + A4 1- has positive, finite and nodegenerate eigenvalues. By the 

expansion 

(116) M(s) = 1 2N + T dC(h) + f][ S dC(h) t dC(t 2 )... f 1 ' 1 dC(t,) 

Jo l=2 Jo Jo Jo 

we see that M is certainly almost surely bounded and that for any < s < 1 

(117) M*_ + (s)M-+(s) = b(s)b(s)* + K{s) 
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with K(s) G C NxN hermitian and E {\\K(s)\\) < Cs 3 ' 2 , while 



(118) b(s)b(s)* = sb(l)b(l)\ 



Now once we see that b(l)b(l)* has almost surely positive and nondegenerate eigen- 
values, we can denote the minimum of eigenvalue spacing and distance to zero of 
spec (b(s)b(s)*) by 5(s), which is O(s), and get 



(119) limP({||AT(s)|| < <5(s)/2» = 1 



by a Markov estimate. Then perturbation theory for the spectrum of hermitian matrices 
(El, Chapter V., Theorem 4.10) proves Lemma |6j 

To see that b(l)b(l)* is positive and nondegenerate almost surely we proceed as in 
the proof of Lemma 2.5.5. in [2J. Let p be the characteristic polynomial of b(l)b(l)*, 
and q(\) = Ap(A). b(l)b(l)* has a degenerate or zero eigenvalue if and only if q has 
a multiple zero, which is when the discriminant D(q) is zero ([2], A.4). D(q) is a 
polynomial in the coefficients of q, and thus a polynomial in the coefficients of b(l), 
which does not vanish identically (as there are of course nondegenerate positive real- 
izations ofb(l)b(l)*). But as the probability density of the entries of b(l) is continuous, 
P (D(q) = 0) = 0, because zero sets of nonvanishing polynomials have Lebesgue mea- 
sure zero. □ 



Thus we know that if we define for any e > the event 
(120) A e = {w G n : T M {e) G D N } 



then P (A e ) — > 1 as e — > 0. For the time being we fix e > and start the equation ( |26| 
on the set A € from s = e with initial data T e (e) = T M {e) to obtain a unique strong 
solution (T e (s; w)) s > e wgj4e . The paths of T e stay inside Djy for all s > e almost surely 
on A € by Theorem |2] To define T e for all times and on the whole probability space, 
let T e (s;ui) = T m (s;uj) if0<s<eorw^ A e . T e is a process with almost surely 
continuous paths. Of course, T € gives rise to processes A| = (T^)" 1 . From the well- 
posedness of the DMPK equation for T £ on A € for s > e, we see by an application of the 
Ito formula that A e solve d25jl on A e for s > e with initial condition A e (e) = X M (e). As a 
consequence, the following N x N matrix-valued processes are well-defined and have 
continuous semimartingale entries: On A e , start with Lf(e) = 0, Rf{e) = 0, and let for 
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s > e 
(121) 

V x l - x ^^ + \l x l - 2A„-1 

dlrjk = da +iAU/ + + tg^ - — ==dlmb 

A " ~ A H / 2y / A^ - X v 



vv 



y/Xl - \ v db Vfl + J XI - X^db^A 2X U -1 

dL = da-^ u + Srfv - i5^—-===d\mb vu 

-v ~ a a« / 2 VK ~ x » 

, p+ , | y/A„(A„ - l)db^ + y%(A„ - gdjj^A 1 r 

d ^ = I rrrf ) ~ / " / 2 v /A^TT 

where we suppressed the e-dependence of X,L,R for a moment. Note that all four 

matrices are skew-hermitian, and furthermore, for (3 = 1, also Lf = and Rf = ■ 
Now decompose A4(e) as in and define matrices U±(s), V±(s) for s > e on A E as 
the unique strong solutions of 

dUUs) = UUs) (dLf(s) - hLt(s)dLt(s)*) , UUe) = U±(e) 

(122) V 7 

dVl(s) = (dRt(s) - 2 dRt(s)*dRt(s)j Vl(s), Vl(e) = V ± (e). 

Existence and uniqueness of these solutions can be easily obtained as long as the de- 
nominators in the coefficients of (121 1 do not get too singular, which in turn can be 
controlled by the stopping-time argument introduced in the proof of Theorem [2] Fur- 
thermore, we have 

Lemma 7. All four matrices U±, V± are unitary. Also, for (3 = 1, U+ = UL and V+ = VI. 

Proof. For the first statement, consider, say U±. U±(e) = U±(e) is unitary by definition. 
On A t , Ito's formula yields for all s > e 

d (UL_UL*) = U%_ (dLf - ^dLfdLf*^ UL* + U e ± (dLf* - ^dLfdLf*^ U£ 

(123) + U^dLfdLf*U^ 
= 

by the skew-symmetry of Lf. In addition, if (3 = 1, we have 

(124) Uj{e) = U^e) = U. (e) = UL (e) , 
and, by lit = L~, 

(125) dUj = U% [dLf - ^dZfdL+*^ = U% (dL; - ^dL~dL~* 
so UL solves the SDE for UL. The proof for the V± matrices is analogous. □ 



22 M. BUTZ 

On A e for s > e we now write 

e^_f A!(g) v^M 

* lSj - V\/A e (*) 2 - 1 A £ (*) 

with A e a diagonal matrix with entries A e uu = i/A^. Now we are ready to replace M by 
the process M e with 

M e (s;uj) = M(s;uj) for < s < e or oj 4 A e 

(127) ~ ~ 

V ; M'^oj) =W(s;uj)S t (s;uj)V e (s;uj) otherwise. 



Also, let 

(128) C(s) 



a + (s) b(s) 
b*(s) a-(s) 



Now as U e (s), S € (s), V e (s) are all well defined on A e as semi-martingales for s > e, we 
can calculate the dynamics of M. e {s) with the Ito formula. In the calculation below, all 
equations hold as stochastic differential equations on the event A e starting from time e. 
Suppressing the e- and s- dependence, denote 

so that 

(130) dU = U [dA + * dAM^ dV = (dB + \dBdB^ V, 

where we used the skew-symmetry of A and B. As dS is a matrix consisting of four 
diagonal blocks, we only give the zv-th diagonal element of each block, reading 

/ dA„ dA„ \ -, / \ v — 1 A„ \ 



After splitting S into its continuous local martingale part Si and a finite variation part 
S2, we have 

< 132 » (^^."(dR*. 
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1 + 



P(N-l) + 2 



2X P X U — X u — X, 



X,, — X, 



ds 



2A„-1 

l 



ds 



/3(JV-l) + 2 







p-£v \ v -\ p 







ds 



/ 



where a factor was included for later convenience. Comparing (132 1 to (121 I, it is 
clear that the bracket processes d*4d5 = dBdS = vanish, and thus the Ito formula for 
M e reads 



(134) 



dM e = d{USV) = (dU)SV + U{dS) V + US(dV) + (dW)S(dV) 
= U(dA)SV + U{dSi)V + US{dB)V 

+ U(dS 2 )V + hu{dAdA)SV + -US(dBdB)V + U(dA)S(dB)V. 



By a straightforward computation using ( |121| and (132 1, we see for the second line 

(135) U{dA)SV + W(dSi)V + US{dB)V = U{dC)SV. 
For the third line of ( [1341 , we have the following trick 

(136) {dA)SdB = (dA) (d<Si + SdB) = (dA) (cLC - d„4)S 
by ( 135^ and analogously 

(137) S(dBdB) = (d«Si+«SdjB)(cLB) = (d£ - d„4)Sd,6 = (dl-dA)(d£-dA)S-dldS 1 . 
So far, 

(138) 

dM e = U(dC)SV 

+ \u{dCdC)SV + U [ {dS 2 )S- 1 + hdAdC - dCdA) - ^(dS^S -1 ) SV. 
By definition of a + , d_ , b, d£d£ = 0, while 



(139) 



- ^(ds^s- 1 



2/3(N- l)+4 V 1 n 







ds 



24 



M. BUTZ 



and 
(140) 

1,. ^ 1 /YdL+ \ /0 db\ / db\ /dL+* 

;(dAd£-d£d.A) = - _ , r _ ~ n ) + ' 



/ n 2A„-1 , \A|E^ 

/3(iV - 1) + 2 2A,-i i ^ V^Fa: n 

\ 2 ^=a: + P2 -p^ A„-A p u 

The matrix in the last line consists of four diagonal blocks again, and we have only 
given the u-th diagonal elements. 

So all terms in ( 138} cancel except for one, and we have on A e , for s > e 

(141) dM € (s) = U e (s)dC{s)S e {s)V e (s) = U e {s)dC(s)U e *{s)M e (s). 
Define the process C e by 

C e (s; u) = C(s; uS) for 0<s<eoru;^^4 e 

(142) f s 

C e {s) = C{e) + J W{t)dC{t)U e *{t) otherwise. 

The stochastic integral in the second line is well-defined by the construction we have 
followed so far, so C e is a well-defined continuous local martingale, and in view of ( 141} 
we have on f2 and for all times s > 

(143) dM e {s) = dC e (s)M e (s). 

To verify that the process A4 e has the same law in the space of continuous functions as 
M., it is enough to show 

(144) (£m>o = (£00W 

For s < e, we even have £ e (s) = C(s), and for s > e, we need to prove that the paths of 

(145) £ e (s) - £ e (e) = l(A e ) J* U e (t)d£(t)U e *(t) + l(A c e )(£(s) - £(e)) 

are distributed as C(s) — C{e), and independent of (£ e (r)) r<e . For the first statement, 
it is enough to notice that the bracket processes of C e and C coincide by the invariance 
property 

(146) W{t) (C(t + 5)- C(tj) U e *{t) = lit + 5) - C{t) 



for all 5 > 0, which follows from the distribution of o+, a_, b. For the independence, it 
is important to note that U e (t) do depend on (£ e (r)) r<e by their initial condition, but 
this dependence factors out, again by ( |146} . So finally we have 

(147) (^m>o = (M*)W 

The process T e is by construction the transmission eigenvalue process of M e , and thus 

(148) (T e ( s )) s > Q i(r^( s )) s>0 , 
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on path space. T M and T e are almost surely continuous processes, as Ai and M e are 
and singular values of a (sub)matrix depend continuously on the matrix. As the "largest 
part" of T e is already defined as a solution of the DMPK equation, we are now in a 
position to compare the distributions of T M and T, and thus to prove Theorem [5] 

Proof. We will prove that the distribution of (T e (s)) s>0 on (C k jv[0, 00), B (C r jv[0, 00))) 
converges to that of (T(s)) s>0 , as e — > 0. As we know that all T e are equal to T M in 
distribution on path space, this can only mean that the distributions of (T M (s)) s>Q and 
(T(s)) s>0 already coincide, which proves Theorem|5| 

From the initial condition for T, its boundedness and its continuity, we directly have 
(149) E[ sup ^ll-TfcOOl^ -+0 (e^O), 



y°^k=l 



0). 



and, recalling that T e (s) = T M (s) for s < e, and thus e-independent, bounded, and 
continuous, also 

(150) E( sup V|l-T fc e ( S )| ) ^0 ( 
Thus we already know 

(151) E ( sup \ T "( S ) " T k( s )\] ^ ( e -> 0). 

For a given e > 0, define the auxiliary process ( f € (s) ) as the unique strong solution 



s>e 



of a DMPK equation driven by B^{s), k = 1, N, s > e, starting from T^(e) = 1 for 
all k, which exists by Theorem [4] On the other hand, for any e > 0, T e (e) G on the 
event A t , so we have the ordering > T e (s) for all s > e on A t by Lemma [3] (here 

and in the following we ignore the singular initial condition for T e , as we have seen in 
the proof of Theorem|4]how to approximate T e from below). Thus on A e , we have 



N N 



(152) Ue (s) = 2 \n( S ) - n(s)\ = y, (n(s) - n( S ) 



k=l k=l 



for all s > 0, and can argue as at the end of the proof of Theorem |4j that on event A e , u e 
is a solution of 

(153) u e (s) = u e (e) + (fm,Tm) dt + [ a T k(t)) 
with a Brownian motion B e and 

f(ffc),Tm) <cu e (t), 

(154) 

with C, C only depending on N. Thus we have again 

(155) E (u e (a) • l(A e )) < E (u e (e) ■ l(A e )) exp(C(s - e)). 
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Fixing an t > e, with an application of the Doob inequality 

E ^ sup (l(A e ) ■ J\{f^r),T^T)) dB e (r) X 

< 4E ( (l(A e ) ■ f a(f^r),T^r)) dB*(r) 



(156) 



< 4C 2 E M(^ e ) • J u £ (r)dr 

< 4C 2 E(l(A e )u e (e)) J exp(C7(r - e))dr 
AC 2 

<^e Ct E(l(A e )u e (e)) 



and the estimate for /, we even find a constant K(t, N) such that 

(157) E ( sup u e (s) ■ l(A e ) ) < K(t, N) (e (1(A £ ) • u e (e)) + y/E (l(A e ) ■ u e (e)) 

\se[e,t] J v 



Together with Lemma [6J and our control of u e (e) from (150 1, 
(158) 

E ( sup V \U( S ) - n(s) ) < K(t, N) (E (u £ (e)) + +N(l-¥ (A e )) -> 

as e — )■ for any fixed t. By a similar argument, (except that we do not have to control 
a possible ill-behaved event Afj, 



(159) 



e( sup Y,\T € k (s)-T k (s)) 







as e — > 0, t fixed. But the combination of ( |151) , (1581, ( 159| l is equivalent to 
(160) E(dist (T e ,T)) -> 

as e — > 0, with dist denoting the metric of local uniform convergence on C r jv [0, oo). A 
fortiori, T e converges to T in distribution on C r n [0, oo). □ 

2.3. General stochastic processes with Coulomb repulsion. In general, we expect the 
following to be true for stochastic processes with Coulomb repulsion. 

Let N G N, and choose bounded qy : 1^ —> M. twice differentiable with bounded 
derivatives, v : R N ^ R N Lipschitz, D k : R N -)• E Lipschitz for all A:,/ = l,...,iV. 
Moreover, let the following symmetry conditions be satisfied: 

C a (k)a(l)( x l, —,%N) = CM —X(j(N)) 
Va(k)(%U -,XN) = V (2^(1), ...X CT (jV)) 

for all permutations a and all a; G R , and assume 



(161) 



(162) 



1 



c«(x) + Q fc (x) > - (L>2(x) + (x)) 
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for all k, I = 1, N, and all x £ Mr with \x k — xi \ sufficiently small. 

Then, for any initial data xi(0) < ... < xn(0), the stochastic differential equation 

(163) dx k = D k (x)dB k + v k (x)dt + Y Cfc;(x) dt 

U Xk ~ xi 

has a unique strong solution (X t ) t>Q , with xi(t) < ... < a?jv(t) for all t > almost surely. 

If, even more, c k i(x) > for all k, I, x, one can also start from singular intial conditions 
a?i(0) < ... < xn(0) and still have a unique continuous stochastic process (X t ) t>0 such 
that for all positive t the particles are strictly ordered x\(t) < ... < xjy(t), and (X t ) t>0 is 
a strong solution to ( 163} . 



A short comment on the assumptions is in order. While ( 162} makes sure that the 



diffusion is controlled by level repulsion at all times, and can essentially not be relaxed 
without allowing for collisions, ( |161} is just one of many ways to make sure that the 
'ODE part" of ( 163} does not allow for intersecting trajectories. We choose to write 



down these assumptions as symmetries, as those should arise quite naturally whenever 
the x k are eigenvalues of a matrix-valued process with invariance properties similar to 



(141. 
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